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NATIONAL ADVISOR! COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 2151 

ESTIMATION OF THE DAMPING IN ROLL OF SUPERSONIC - 
LEADING-EDGE WING-BODY COMBINATIONS 
By Warren A. Tucker and Robert 0. Piland 

SUMMARY 


A description is given of a linear-theory method for the estimation 
of the damping in roll of supersonic-leading-edge wings mounted on cylin- 
drical bodies. The error resulting from an approximation essential to 
the method is shown to be small. 

The damping in roll is calculated for configurations having rectan- 
gular and triangular wings) for most of the range of variables the net 
effect of the body is to decrease the damping in roll. The body effect 
is small for body diameters less than about 30 percent of the wing span 
but increases rapidly with body size above this value. 


INTRODUCTION 


The damping in roll of wings of various plan forms in supersonic, 
flight has been the subject of severed, investigations. The theory for 
the isolated wing is well-advanced) in connection with the present work 
the results for the rectangular wing (references 1 and 2) and the trian- 
gular wing (reference 3) are of interest. In actual aircraft, however, 
the wing is commonly mounted on a body so that a method of evaluating 
the effect of the body is desirable for the purpose of applying isolated- 
wing theory to practical cases. 

The present paper describes a linear-theory method for estimating 
the damping in roll of thin super sonic -leading-edge wings mounted on 
cylindrical bodies. The method is applied to the calculation of the 
damping in roll of configurations having rectangular and triangular 
wings. The triangular-wing configuration with a subsonic leading edge 
has been treated by other investigators in reference 

The accuracy of the method in determining the pressure distribution 
in the region of the wing affected by the presence of the body is diffi- 
cult to assess, but within the range of applicability of the method 
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(small wing thicknesses) no viscous effects) the error in total rolling 
moment can .he ‘shown to he small. 


SYMBOLS 


s 

a 


h = 2s 
h 


M 

p = 

ph/2V 

S 

L 


distance from body center line to wing tip 
radius of cylindrical body 

rolling angular velocity, positive when clockwise viewed 
from rear 


distance from body center line to limit of juncture region 
(see fig. l) 

free-stream velocity 

free-stream Mach number 

1 

wing -tip helix angle 
free-stream dynamic pressure 

wing area, including the part of the wing covered by the body 


rolling moment, positive when acting clockwise viewed from 
rear 

rolling-moment coefficient 

'Be, 


damping-in-roll derivative 


27/ 


Kl 


A 

*,y 


damping-in-roll derivative for isolated wing 
aspect ratio (b^/s) 

coordinates for field points (see fig. 2) 
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M 

G 

0 

d ^0 

0x - as 


coordinates for source points (see fig. 2) 
local source strength 
disturbance -velocity potential 


1 

m = k 


slope of leading edge (see fig. 2) 


pressure difference between upper and. lower surfaces 


pressure coefficient 


ASALTSIS 


The basis of the analysis is the representation of the rolling wing 
and body by a distribution of sources chosen so as to satisfy approxi- 
mately the boundary conditions of the problem. After this distribution 
has been chosen, the determination of the resulting pressure distribu- 
tion and rolling moment is straightforward, although somewhat lengthy. 

Source distribution .- The geometry of the rolling wing-body system 
is shown in figure 1(a). The part of the fuselage to which the wing is 
attached is assumed to be a cylinder of radius a, and the whole system 
is rolling at an angular velocity p. 

The right wing, for example, can be replaced by a distribution of 
point sources the strengths of which are equal to l/:rt times the local 
vertical velocity (reference 5)- For the present case, then, the source 
strength is directly proportional to the distance from the body center 
line (see fig. 1(b)). 

The boundary condition to be satisfied on the body is, of course, 
that there be no flow normal to the body. In order to satisfy this con- 
dition, the source distribution necessary to represent the body must be 
such that it acts together with the source distribution for the wing to 
produce a flow which has no component normal to the body. The determi- 
nation of such a source distribution would be a very difficult task, and 
the distribution when found would very likely be so complicated as to 
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present further difficulties in the succeeding steps. An approximation 
is therefore in order. 

Consider the case for which the body radius is only slightly less 
than the distance h (see fig. 1(a)). For this condition the effect of 
the body on the wing will be much the same as that of a large flat plate 
mounted normal to the wing along the wing-body juncture. The source 
distribution req uir ed to represent this flat plate is simply that shown 
in figure l(c). As the size of the body relative to the distance h is 
decreased the effect of the body becomes less and less like that of a 
flat plate until finally the source distribution is simply that for the 
wing alone (see fig. 1(d)). 

This knowledge of the conditions at the end points has been used 
to select a source distribution which represents the effect of the body 
in an appro ximat e manner. The chosen distribution is shown in fig- 
ure 1(e) j the slope of the distribution curve varies linearly with a/h. 
This distribution is at best a crude approximation to the actual state 
so far as the bo undar y condition on the body is concerned. However, as 
shown in the section entitled "Discussion, " the resulting error in 
rolling moment is small. 

Pre ssure di stribution . - The source distribution being chosen, the 
method of reference 5 can be used to find the pressure distribution in 
the wing area between the wing-body juncture and the Mach line from the 
leading edge of the juncture. This area is hereinafter called the junc- 
ture region. This region is the only one in which the pressure need be 
determined and in which the body effect 1 b felt directlyj the pressures 
acting on the body are quite ineffective in producing rolling moments 
and the pressures over the other areas of the wing are known from refer- 
ences 1 to 3 . 

The source distribution -of figure 1(e) is sufficient to show that, 
at least for the rectangular -wing case, the pressure in the juncture 
region has a greater absolute value than would be the case were the body 
not present. The corresponding increase in damping, however, is opposed 
by the decrease which results from covering part of the wing with the 
body. The relative magnitudes of the two effects cannot be evaluated 
from qualitative reaso nin g. 

The axis system for the determination of the pressure in the junc- 
ture region is shown in figure 2. The disturbance -velocity potential 
at a point (x,y) is given by ' 


0 = - fT Gf(|,ii) cLtj 


JJ 


( 1 ) 
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where G is the local source strength, the integration is over the 
entire shaded area in figure 2, and 




(* - 


£) 2 - 


P 2 (y - n ) 2 


-1/2 


The integration is most conveniently carried out over the three sfeparate 
shaded areas in figure 2 so that equation (l) becomes 


* f = 

P r 


pgy+x 

p+k 


v>c+py-pq 


Jy 


py Px 
JO Jkt] 
ro 

Py-x 


I hr] 

•py+pT) 


(a + Tj)f(£,ii) dS drj + 


(a + Tj)f( |,ti) dS dn + 


''x-py+pTj 

/ 

0 a ->\~1 

-*n 

! " ( 

2 h - Vll 


( 2 ) 


The first integration with respect to | gives 



PfZt* 

p+k 


uo 


(a + ti) 


cosh 


-1 


X - fa) 

P(y - tJ) 


dt] + 


’0 




[• 


p+k 




cosh 


-1 


x + kr| 

p(y - n) 


dii 


( 3 ) 


The partial derivative is necessary for the determination of the 

pressure. [This derivative is found most conveniently at this point in 
the analysis by differentiation under the integral sign. . This operation 
results in 
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i 0x = 


' py+x 

p+k 


(a + q) dg . 


0 Vx 2 - p 2 ;? 2 + (2p 2 y - 2kx)T) - (p 2 - k 2 )q 2 

> [f - ( 2 § ~ X ) j] d 4 

py-x \/x 2 - p 2 y2 + (2p2y + 2kx)Tj - (p 2 - k2)q 2 
p+k 




These integrals are easily evaluated. If a pressure coefficient is 
defined as 



(5) 


where P is the pressure difference between the upper and lower sur- 
faces^ then 

Cp -= ft* (6) 

Substitution of equation (4) into equation (6) results in 



'/p 2 - k 2 


, — -1 kx - p 2 y 
cos p(7-r^y + 


cos 


_1 kx + 


fC 


X + 


ky)| 


P^ - kx 

( s 2 - ^) 3/2 


cos 


-1 kx - P 2 y 
P(x - ky) 




P^y + kx 
0 2 - k2) 3/2 


COS 


-1 kx + p^ 
P (x + ky) 


(T) 


This equation gives the pressure coefficient acting over the region of 
the wing affected by the bodyj as has been mentioned before, the pres- 
sures acting over the other areas are known from references 1 to 3. 
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For the rectangular -wing case, k is put equal to zero. The 
quantity ^2 ^ - lj can be expressed in terms of more familiar parame- 
ters. For the rectangular -wing: 

a 2_ 

o a 1 _ s ~ pA 

h a _2_ 

s + PA 

and for the triangular wing: 



Damping -in-roll derivative .- The damping -in-roll derivative C2 p 

is found for either the triangular or the rectangular plan form by suit- 
able integration of the various pressure distributions. The rather, 
lengthy integrations are not included here. The wing area and span upon 
which the coefficient is based include those parts of the wing masked by 
the body. 

The result for the rectangular -wing configuration is given by the 
following equation: 



_ 2 1 1 1 

3 PA 3 (PA) 2 12(PA) 3 


2 - 
£ s 

, 2/a' 

£ la 

k 

(M) 

W PA\s t 

' ” (PA) 2 8 " 

3*(PA)3_ 


(8) 


This equation reduces correctly to the equation for the isolated wing 
(references 1 and 2) if a/s is put equal to zero. The equation is 
valid whether or not the Mach lines cross each other, so long as the 
Mach line from the leading edge of the wing-body juncture does not cross 
the wing tip. This condition can be expressed mathematically by the 
expression 



(9) 
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f T T ViR corresponding result for the triangular -wing configuration is 
given by the following equation: 




Substitution of — = 0 causes all but one term to vanish. The value l/3 
. s 

given by this term is the correct value for the isolated wing (refer- 
ence 3). The only restriction on equation (10 ) is that the Mach line 
from the leading edge of the wing-body Juncture must be behind the 

leading edge of the wing (fij~ = As PA approaches infinity, equa- 

tion (10) becomes 



( 11 ) 


DISCUSSION 


Effect of approximate source distribution .- In the preceding 
analysis, the presence of the body has been represented by a rather 
arbitrarily selected source distribution. The error arising from this 
approximation should therefore be examined so that the validity of the 
results can be assessed. One method of evaluating the error is to 
determine how closely the boundary condition of no flow normal to the 
body is met by the velocity field induced by the source distribution of 
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figure 1(e). If this approach is adopted, it is not difficult to show 
that the maximum value of Vq/V is of order (pb/2V)(a/s), where vq 
is the velocity normal to the body. However interesting this result 
may be from an academic point of view, it supplies no direct answer to 
the more practical question of how much error can be expected in the 
estimated value of C 2 p for a given wing on a given body. 

The maximum possible error in C2 p can be determined by the con- 
sideration of two extreme cases. As one extreme, the slope of the source 
distribution representing the effect of the body can be taken equal 
to p/rt, regardless of the body size (see fig. 1(e)). This case will 
give the maximum possible value of C 2 p (the physical interpretation 

is that the effect of the body is taken as that of a large flat plate 
mounted normal, to the wing along the wing -body juncture). As the other 
extreme, the slope of the source distribution representing the effect 
of the body can be taken equal to - p/it, again without regard to body 
size. This case, which will result in the minimum value of C 2 p , can 

be regarded for the rectangular-wing configuration as an assunption that 
the body has no effect on the part of the wing outside the body. As 
simple a physical interpretation is not possible for the triangular -wing 
case. The difference between the two values of C 2 p will of course be 

the maximum possible error. 

The pressure distribution in the juncture region for these two 
extreme cases is shown in figure 3 for a particular set of rectangular 
wing -body combinations, together with the pressure distribution given 
by the present approximation. The inference is clear that large errors 
in span load distribution are possible, particularly at stations close 
to the wing -body juncture. The error in C 2 p , however, is not neces- 
sarily large, because of the relative nearness of the juncture region 
(compared with the outer parts of the wing) to the moment axis (body 
center line). That the error in C 2 p f° r "the rectangular -wing case is, 
in fact, acceptably small for engineering purposes is demonstrated by 
the curves of figure 4 , which is self -explanatory. The information cor- 
responding to that shown in figure 4 for the rectangular-wing case is 
given in figure 5 for the triangular wing} the error in C2 p is again 
small. 

Rectangular wing .- Now that the error introduced by the use of an 
approximate source distribution has been shown to be small, the results 
of calculations made with the use of equations (8), (10), and (ll) may 
be discussed. The effect -of body size on the damping in roll of 
rectangular -wing configurations is shown in figures 6 and 7* The ques- 
tion posed in the section entitled "Pressure distribution," namely 
whether the effect of covering part of the wing with a body or the 
effect of the increased damping in the juncture region is pred ominant . 
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is here resolved) the effect of covering part of the wing is the greater 
for most of the range of variables. The net effect of the body is not 
very important, however, for values of a/s less t han about 0.3« The 
curves of figure 7 show that the percentile effect of the body is not 
greatly dependent upon PA. 

Triangular wing .- The effect of body size on the damping in roll of 
triangular-wing configurations is shown in figure 8. As was the case 
for the rectangular wing, the net effect of the body is generally to 
decrease the damping in roll, although some opposite effect is present 
for small values of a/s. 


The percentile effect of the body for the triangular -wing case is 
shown in figure 9* which corresponds to figure 7 for the rectangular- 
wing case. In contrast to the rectangular -wing case, the value of PA 
has an appreciable influence on the magnitude of the body effect. The 


curve for ^ — >0 is taken from reference. % the value — *0 

refers physically to the condition that the wing lea din g edge is highly 
swept back relative to the Mach lines. 


Other considerations.- Although the damping-in-roll derivative C^p 
has been worked out only for rectangular and triangular plan forms, 
the pressure coefficient in the .Juncture region given by equation (7) is 
equally applicable to any other plan form with a straight supersonic 
leading edge. There is also little doubt that the results may be 
applied, within the accuracy of the present analysis , to configurations 
having bodies which are not true cylinders, but which are almost cylin- 
drical (quasi -cylindrical). 


There are at least three limitations Inherent in the analysis. The 
first is, of course, the use of an approximate source distribution to 
represent the body effect) the attendant errors have been discussed in 
the first part of this section. The second is the assumption of van- 
ishingly small wing thickness. The effect of moderate thickness in the 
absence of viscosity may be expected to be small (see, for example, 
reference 6). The third, and possibly most severe, limitation is the 
assumption of an inviscid fluid. The magnitude of this error is practi- 
cally impossible to evaluate theoretically, and there is not yet much 
applicable experimental data. 


CONCLUDING REMARKS 


A linear-theory method for the estimation of the effect of a 
cylindrical body on the damping in roll of wings with supersonic lea din g 
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edges is presented. The essential feature of the method is the repre- 
sentation of the effect of the body on the wing by a rather arbitrarily 
chosen source distribution. The error in damping in roll ■which could 
result from this procedure is shown to be small. 

The dangling in roll is calculated for configurations having 
rectangular and triangular wings. The body effect is found to be small 
for body diameters less than about 30 percent of the wing span but 
increases rapidly with body size above this value. 

The results for the supersonic-leading-edge triangular- wing case 
are compared with those obtained by other investigators for the subsonic- 
leading-edge case. The leading-edge condition is shown to be an 
important parameter. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., April 12, 1950 
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(e) Final source distribution. 
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Figure 1.- Steps in choice of source distribution. 
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Figure 3 -- Various 


approximations to pressure in juncture region 


- = 0.3; PA = k. 
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‘ Rectangular wing-body combination. 







Figure 4.- Effect of various approximations to pressure in juncture 
region. Rectangular wing-tody combination. 
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Figure 9.- Relative effect of "body for triangular-wing conf igurations . 
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